A new class of oscillatory radial basis functions

Radial basis functions (RBFs) form a primary tool for multivariate interpolation, and they are also receiving increased attention for solving PDEs on irregular domains. Traditionally, only non-oscillatory radial functions have been considered. We find here that a certain class of oscillatory radial functions (including Gaussians as their limiting case) leads to non-singular interpolants with intriguing features especially in the different limit in which they are made increasingly flat. This latter limit is important in that it generalizes traditional spectral methods to completely general node layouts. Interpolants based on the new radial functions appear immune to many or possibly all cases of divergence that in this limit can arise with other standard types of radial functions (such as multiquadrics and inverse multiquadratics).